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Abstract 

The CHORUS experiment, designed to search for oscillations, consists of a nuclear emul- 

sion target and electronic detectors. In this paper, results on the production of charged particles in 
a small sample of charged-current neutrino- and anti-neutrino-nucleus interactions at high energy 
are presented. For each event, the emission angle and the ionization features of the charged particles 
produced in the interaction are recorded, while the standard kinematic variables are reconstructed 
using the electronic detectors. The average multiplicities for charged tracks, the pseudo-rapidity dis- 
tributions, the dispersion in the multiphcity of charged particles and the KNO scaling are studied in 
different kinematical regions. A study of quasi-elastic topologies performed for the first time in nu- 
clear emulsions is also reported. The results are presented in a form suitable for use in the validation 
of Monte Carlo generators of neutrino-nucleus interactions. 
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1 Introduction 

In the study of multi-particle production processes, the multiplicity of charged particles is an im- 
portant global parameter reflecting the dynamics of the interaction. Different phenomenological and the- 
oretical models can be tested and therefore multiplicities have been studied in experiments with different 
particle beams, making use of various techniques and over a wide range of kinematic variables. 

The characteristics of charged-particle multiplicity distributions have been studied in detail in 
high energy hadronic collisions [IJ, in e~e~^ annihilation EllSl, and in interactions of neutrinos and 
anti-neutrinos on nucleons and light nuclei ||4l|5l|6l|7l[8l. In particular, Bubble Chamber experiments I9l 
[TOl [Tn performed measurement of the charged particle multiplicity, produced in (I^/i)n and 0^^)p 
interactions. They also studied the multiplicity moments as function of kinematical quantities. 

The sub-micron spatial resolution of nuclear emulsion allows both the investigation of the event 
topology and the measurement of the angular distribution of charged particles to be performed with high 
accuracy. Therefore, it is well suited to study particle-production multiplicities. However, such data on 
neutrino interactions in nuclear emulsion are relatively scarce llT2l[T3l and don't exist for anti-neutrino 
interactions in nuclear emulsion. 

In this paper, we report on a study of charged-particle multiplicities produced in high-energy 
charged-current (anti-)neutrino interactions in a nuclear emulsion target measured with the CHORUS 
hybrid detector. Such data are useful in tuning interaction models in Monte Carlo event generators. We 
also investigate the KNO scaling lfT4]| behaviour in different kinematical regions. 

2 The experimental apparatus 

The CHORUS detector ifTSi is a hybrid setup that combines a nuclear emulsion target with several 
electronic detectors: trigger hodoscopes, a scintillating fibre tracker system, a hadron spectrometer, a 
calorimeter and a muon spectrometer. 

The target consists of Fuji ET-7B |[T6l nuclear emulsions with properties given in Table [T] The 
nominal sensitivity is 30 grains per 100 ^m for minimum-ionizing particles. The target is segmented into 
four stacks and has an overall mass of 770 kg. Each of the stacks consists of eight modules of 36 plates of 
size 36 cm x 72 cm. Each plate has a 90 ^m plastic support coated on both sides with a 350 ^m emulsion 
layer. Particle track directions are nearly perpendicular to the emulsion sheets. This configuration permits 
fast automatic scanning of the emulsion sheets. 

Each target stack is followed by three interface emulsion sheets with 90 fim emulsion layers on 
both sides of an 800 fim thick plastic base and by a set of scintillating-fibre tracker planes U7] . The 
interface sheets and the fibre trackers provide accurate particle trajectory measurements which are ex- 
trapolated into the emulsion stack in order to locate the neutrino interaction vertices. The accuracy of the 
fibre tracker prediction is about 150 /Lim in the track position and 2 mrad in angle. 

The electronic detectors downstream of the emulsion target include a hadron spectrometer that 
measures the bending of charged particles in an air-core hexagonal magnet ifTSl . a calorimeter where the 
energy and direction of showers are measured and a muon spectrometer which measures the charge and 
the momentum of muons. The energy resolution for hadronic showers of the calorimeter is a{E)/E = 
(0.323 ± 0.024)/y/^/GeV + (0.014 ± 0.007) [ 19 1. The muon spectrometer consists of six magnetized 
iron toroids sandwiched between seven tracking stations (drift chambers and streamer tubes). The muon 
momentum resolution varies from k,\5% in the 12-28 GeV/c range |20| to 19% at about 70 GeV/c |[T5l . 
as measured with test-beam muons. 

The West Area Neutrino Facility (WANE) of the CERN SPS provided an intense beam of neutrinos 
with an average energy of 27 GeV. The beam consisted mainly of with a contamination of ~5% 
and ~1% i/g- The CHORUS detector was exposed to the wide band neutrino beam (section 5 of ifTSl for 
details) of the CERN SPS during the years 1994-1997, with an integrated flux of 5.06 x 10^^ protons on 
target. In this four year exposure more than 10^ neutrino interactions were accumulated in the emulsion 
target. 

3 Analysis 

The emulsion scanning has been performed by fully automatic microscopes each equipped with a 
CCD camera and a fast custom built processor (originally called 'Track Selector') 11211 capable of iden- 



Table 1: Atomic composition and main features of the nuclear emulsions (Fuji ET-B7) used in the CHO- 
RUS experiment IMl 



Element Atomic number 



Mass (%) 



Mole fraction (%) 



Iodine (I) 53 

Silver (Ag) 47 

Bromine (Br) 35 

Sulphur (S) 16 

Oxygen (O) 8 

Nitrogen (N) 7 

Carbon (C) 6 

Hydrogen (H) 1 



0.3 

45.5 

33.4 

0.2 

6.8 

3.1 

9.3 

1.5 



0.06 

11.2 

11.1 

0.2 

11.3 

5.9 

20.6 

40.0 



Mean number of nucleons 
Density 

Radiation length 

Nuclear interaction mean free path 
Concentration of AgBr 



36 protons, 45 neutrons 
3.73g/cm3 
2.94 cm 
38 cm 

45.5% in volume 



tifying tracks inside the emulsions, and measuring on-line their parameters. The track finding efficiency 
of the Track Selector is higher than 98% for tracks with an angle of less than 400 mrad 1,22,1 with respect 
to the beam direction, perpendicular to the emulsion plates. 

The location of the emulsion plate containing the interaction vertex is achieved with a procedure 
which is called scan-back [23l. As a first step, tracks reconstructed by the electronic detectors are fol- 
lowed upstream into the emulsion stacks. A track, extrapolated from the fibre-trackers and found in the 
interface emulsion sheets, is followed upstream in the target emulsion stack using track segments recon- 
structed in the most upstream 100 of each plate. If the track is not found in two subsequent plates, 
the first plate with no track hit is defined as the plate containing the vertex. In total, about 150 000 CC 
events with at least one reconstructed muon in spectrometer, have been located in emulsion as a result of 
this procedure. 

For the present measurement, from the 150000 located events with at least one reconstructed muon 
in the spectrometer, a sample of 1208 events was randomly selected to be visually inspected and to 
measure the parameter associated with the charged tracks at the (anti-)neutrino vertex. We use the sign 
of the charge of the highest-energy muon in the events, as determined by the muon spectrometer, to 
determine whether the interaction was induced by a neutrino or an anti-neutrino. In order to have a larger 
anti-neutrino sample, events located in various modules of the emulsion stacks have been collected, while 
the neutrino sample consists of events located in one half module of an emulsion stack. In 627 (581) of 
the selected events the charge of the leading muon is negative (positive). 

To evaluate the number of genuine z7^ charged-current (CC) interactions, one has to subtract from 
this sample the contamination due to: 

- CC events with the fi^ reconstructed as a /i"*"; 

- punch-through hadrons traversing the 5.2 interaction lengths of the calorimeter and reconstructed 
as positive muons in the spectrometer. This can happen in CC events where the fi is not identified 
or in NC interactions; 

In order to reduce this background, a set of selection criteria improving the reliability of the muon 
reconstruction is applied to the events with a muon reconstructed with positive charge: 

- low momentum muons stopping inside the spectrometer, for which the reconstruction is based on 
only a few hits in the electronic detectors, are discarded; 

- the impact point on the spectrometer entrance surface is required to be within a 150 cm radius 
centered around its axis to ensure a long potential path through the toroidal magnets; 

- the track length is required to be greater than three spectrometer gaps (two magnets) and consis- 
tency checks are made among the various reconstruction algorithms; 
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Figure 1 : Pseudo-rapidity distributions for tracks classified as sliower (full histogram) and grey (dashed 
histogram) by the scanner. 

- the trajectory of the muon measured in the spectrometer is required to match with a track detected 
in the fibre trackers, and a selection criterion is applied on the basis of of the global track fit. 
These criteria select 529 events with a reconstructed positive muon. The contribution from vr and 
K decays into a ^, reconstructed in the spectrometer is negligible, as well as z/^ induced charm events 
with the reconstructed as a The remaining background to the sample of /i+ due to events is 
composed of punch-through tracks of pions and kaons and amounts to 5.5 ± 1.9 (stat) it 0.6 (syst) events. 

We further require that square of the invariant mass of the hadronic system, W"^, of (anti-)neutrino 
events is greater than 1 GeV^/c^. After this selection the number of u^-Ps. (p^-A) interactions is 496 
(369). 

4 Measurement procedure and charged prong classification in nuclear emulsions 

When a charged particle passes through nuclear emulsion it forms a latent image by ionizing the 
AgBr crystals along its trajectory. After a suitable chemical processing of the emulsions such ionization 
sites induce the deposition of Ag atoms as dark spots (grains) thus making the trajectory of the charged 
particles visible. Tracks are usually classified according to their grain density (number of Ag grains 
per unit length) as shower, grey and black prongs following the procedure described in Refs. |[25l l26l . 
Shower tracks correspond to relativistic (or nearly relativistic) singly charged particles with a velocity 
j3 > 0.7. Their grain density is g < 1.4 ( 9o represents the plateau grain density of a highly relativistic 
singly charged particle, which equals 29.6 grains per 100 fim for the CHORUS emulsion). Grey tracks 
correspond to charged particles with velocity 0.25 < /3 < 0.7. They are interpreted as recoil nucleons 
emitted during the nuclear cascade. Their grain density is 1.4 go < g < 10 go. The black prongs 
correspond to particles with a velocity /3 < 0.25. These are produced by low energy fragments (protons, 
deuterons, alphas and heavier fragments) emitted from the excited target nucleus. The associated grain 
density is g > W go- 

The above criteria commonly used in nuclear emulsion experiments are very difficult to apply 
in the CHORUS analysis since the emulsion sheets were exposed perpendicular to the beam direction. 
Furthermore, electronic measurement of the energy and momentum for every charged particle was not 
possible due to geometrical acceptances. 

The following procedure was used for the CHORUS data. Once the vertex plate had been located 
by the automatic system, the event was checked and measured extensively with a manually controlled 
microscope system. The black prongs at a neutrino interaction vertex, due to heavily ionizing particles, 
have short path lengths and usually stop within one emulsion plate so that they can be recognized. For 
the remaining mixture of shower and grey prongs, we measured the particle directions. These particles 
are mainly pions with a small contamination of protons. 



Ambiguity in classifying shower and grey prongs can arise either due to the difficulty in strictly 
applying the ionization criteria in the manual scanning or due to variations in the quality of the optics of 
the microscope. To overcome these limitations we decided to classify these prongs by using the pseudo- 
rapidity variable 1,27 J : 

e 

1] = — in tan — (1) 

where 9 is the emission angle of the prong with respect to the neutrino direction. This has the advan- 
tage of being independent of the scanner and of the microscope optics, allowing us to compare in a 
straightforward manner the multiplicity measurements with the theoretical models. The pseudo-rapidity 
distributions for tracks classified as shower and grey by the scanner, both for neutrino and anti-neutrino 
interactions are shown in Fig. [T] These plots show that a scanner-independent classification, based on 
the pseudo-rapidity, is possible and it is consistent with the traditional one. In the following, all prongs 
with > 1 are classified as shower particles. By applying this cut the fraction of shower (grey) tracks 
miss-identified as grey (shower) is (1.23±0.19)% ((7.27±1.40)%). The multiplicities of shower, grey 
and heavy (grey-i-black) prongs are denoted by ng, Ug and n^, respectively. The total number of charged 
hadrons classified as shower particles in an event is defined to be rich = "-s — 1> namely the number of 
shower tracks minus the muon track. 

5 Efficiency Evaluation 

Reconstruction and location efficiencies were evaluated making a detailed simulation of the detec- 
tor response using a program based on GEANT3 [28 1. Large samples of deep-inelastic neutrino interac- 
tions (DIS) were generated using the beam spectrum from the CHORUS JETTA generator ll29l based on 
LEPTO |[30l and JETSET OTl . Quasi-elastic (QE) interactions and resonance production (RES) events 
were generated with the RESQUE |[32]| package with a rate of 9.6% relative to deep-inelastic scattering 
reactions in the neutrino case (26% in the anti-neutrino case). The simulated response of the electronic 
detectors was processed through the same reconstruction program as that used for the experimental data. 
The location efficiency was parametrized by a function of the primary muon momentum and angle on 
which it depends only weakly. The reconstruction and location efficiency as a function of the invariant 
mass of the hadronic system and of the shower prong multiplicities rich is given in Tables [2] and [3] 
for v-A and u-A interactions, respectively. The samples are normalized in such a way that the efficiency 
at the first bin is taken as 1.00. Given the fact that Bjorken y distributions in neutrino and anti-neutrino 
interactions are different, the momentum distribution of the positive muon from CC interactions is 
harder than that of the negative muon in CC interactions. Therefore, the reconstruction and location 
efficiencies of CC events are higher than that of CC events. The efficiency does not depend on the 
black and grey track multiplicity at the primary vertex. 

In (p^) CC interactions, the full kinematics of the event can be reconstructed from the measure- 
ment of the muon momentum p^, the angle 6^ of the produced muon with respect to the beam axis, and 
Efiad, the energy transfer to the hadronic system: 

Eu = Efj_ + Ehad 

QI = 2Ey{E^ - Pf, cosO^) - 
Wl = 2mn{E, - E^) + ml-Ql 

where Ey and E"^ are the energy of the incoming (anti-)neutrino and (anti-)muon energy respectively, -Ql 
the squared four-momentum transfer, m,„ and are the mass of the nucleon and muon respectively, W"^ 
is the square of the invariant mass of the hadronic system. The average neutrino energy of the neutrino 
sample is {E^) = (38.0 ± 1.4) GeV, the mean-square momentum transferred to the hadronic system 
is {QD = (8.6 ±0.5) (GeV/c)^, and the mean-square of the invariant mass of the hadronic system is 
= (26.2 ± 1.3) GeV^/c^. For the antineutrino sample these values are {Ejj) = (42.5 ± 1.6) GeV, 
(Q^) = (4.6 ± 0.3) (GeV/c)2, and = (17.7 ± 0.8) GeVVc'', respectively. All distributions presented 
hereafter include the efficiency correction. 



6 Multiplicity Distributions 

The hadronic shower and heavy prong multiphcity distributions based on the pseudo-rapidity 
selection are shown in Fig. |2] and the distribution as a function of emission angle of shower tracks 
is given in Table I^Pl The average number of shower prongs in neutrino-nucleus interactions is 
(nch(i^-A)) = 3.4ib0.1. The average heavy prong multiplicity is measured to be (nh(i^-A)) = 4.7±0.2. 
The average number of shower and heavy prongs in anti-neutrino induced events are measured to be 
(nch(i^-A)) = 2.8 ± 0.1 and (nh(^'-A)) = 3.5 ± 0.2, respectively. These two quantities were never 
measured yet in nuclear emulsion. 

The average shower prong multiplicities rich for v-A and V-A interactions are plotted as a func- 
tions of W"^ in Fig. [3] The numerical values are given in Tables [5] and [6l The mean multiplicities are in 
good agreement with a linear dependence on In W'^, 

{rich) = a + blnW'^ . (2) 

This expression was fitted to the average multiplicities for between 1 GeV'^/c^ and 200 GeY^/c^. 
The values for the fit parameters a and b were found to be 

nch(z^-A) = (0.45 ± 0.24) + (0.94 ± 0.08) In W'^ 

nch(l^-A) = (0.53 ± 0.20) + (0.82 ± 0.08) In W"^. 

The fitted values of a and b obtained in other experiments are shown in Table|7] Our result on a disagrees 
with the unpublished result of (IT] at two sigma level. However we note that the value of a found in 
Ref L12J is much larger than the corresponding value found in all other experiments. 

For the first time a sample of (anti-)neutrino events measured in emulsions is large enough to 
study (QE-i-RES)-like topologies. In order to have a minimum bias sample of (QE-i-RES)-like events, 
the > 1 GeV^/c'^ cut was not applied to the located events. Hence, the starting sample of f^-A 
(i7^j-A) interactions becomes 627 (529). An event is defined as being (QE-i-RES)-like, if the number of 
shower prongs is zero or one and the number grey prongs zero or one for for u-A interactions regardless 
of the number of black tracks. In order to obtain (QE-i-RES)-like enriched sample in V-A interactions, 
the sum of shower prongs and grey prongs is required to be one or zero regardless of the number of 
black tracks. The efficiency of this selection is estimated from the MC simulation to be 72% (81%) 
for u^-A (pi^-A) interactions. In order to suppress background from DIS interactions we further require 
W"^ < 10 GeV^/c'^ for both samples. After this selection the number of (QE-i-RES)-like events is reported 
in Tabled] The average values of kinematical parameters, (E^), {Ql), {W^), corresponding to these 
events are given in Table [sIF^ The number of background events that mimic (QE-i-RES)-like topology is 
obtained from the MC simulation to be 83.4 and 73.5 events for z/^-A and V^-A interactions respectively. 
The ratio of reconstruction and location efficiency of (QE-i-RES)-like events to that of all CC events is 
found to be 1.22 (1.13) for v^^-A (U^-A) interactions. 

After applying the efficiency and background corrections, the fraction of (QE-i-RES)-like events 
is found to be (13.4ibl.0±2.0)% for and (26.3±1.4ib3.9)% for 17^ interactions, respectively. These 
values are in reasonable agreement with the ones used in the Monte Carlo simulation (Section 5). We have 
accounted for a systematic uncertainty of 15% coming from the efficiency and background estimation by 
Monte Carlo modeling. 

The sub-sample of (QE-i-RES)-like (u^) events with neither black nor grey prongs is important 
for the understanding of nuclear mechanisms involving hadrons in the nucleus. The fraction of this type of 
topology is measured as (1.2ib 0.4ib0.2)% for and (9.5±I.0±1.4)% for interactions, respectively. 
For the neutrino events the proton in the final state is absorbed by the nucleus without any visible activity. 
On the other hand, this fraction is larger for anti-neutrinos since in the final state there is a neutron. 



The corrected number of events is 1231(844) for Vf^-A (V^-A) interactions 

Due to detector resolutions, there are 62(52) events with <0 in neutrino(anti-neutrino) sample. For these events, W'^ is 
taken to be . 



Table 2: The reconstruction and location efficiency in u-A as a function of the invariant mass of the hadronic system W"^ and the number of shower tracks. 
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Table 3: The reconstruction and location efficiency in u-A as a function of the invariant mass of the hadronic system W"^ and the number of shower tracks. 
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Figure 2: Multiplicity distributions for hadronic shower (top) and black plus grey (bottom) tracks in 
i^ju-nucleus and 17^-nucleus interactions 

7 Dispersion and Koba-Nielsen-Olesen distributions 

We have also investigated the dependence of the dispersion, Dch, of the shower hadron multiplicity 
distribution, on the average multiplicity. The dispersion D^h is defined as 

i^ch = ^{nl) - . (3) 

For independent particle production, the multiplicity follows a Poisson distribution. Therefore, one gets 

-Dch = VKh). (4) 



Table 4: The number of shower hadrons in (u^) as a function of emission angle 6 where Ntr is the 
number of hadron tracks is the ratio of the number hadrons tracks to the total number of events. 



(Radian) (6) Ntr 



0.000 -1-0.050 0.032±0.001 (0.031±0.001) 346(161) 0.28(0.19) 

0.050 H- 0.100 0.075±0.001 (0.073±0.001) 554(282) 0.45(0.34) 

0.100 -=-0.150 0.125±0.001 (0.123±0.001) 510(270) 0.41(0.32) 

0.150^0.200 0.174±0.001 (0.176±0.001) 394(280) 0.32(0.33) 

0.200 H- 0.300 0.245±0.002 (0.246±0.002) 774(394) 0.63(0.47) 

0.300^0.400 0.350±0.002 (0.347±0.002) 538 (330) 0.44(0.39) 

0.400^0.500 0.45 1 ±0.002 (0.449±0.003) 380(278) 0.31(0.33) 

0.500^0.600 0.544±0.003 (0.544±0.003) 290(171) 0.24(0.20) 

> 0.600 0.65 1 ±0.003 (0.655±0.004) 210(121) 0.17(0.15) 



Table 5: The shower multiplicity distributions in v-K as a function of W"^ (errors shown are statistical only). 
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Table 6: The shower multipUcity distributions in v-A as a function of W"^ (errors shown are statistical only). 
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Figure 3: The hadronic shower prong multiplicity distributions as a function of In W for v-A, and u-A 
interactions. 
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Figure 4: The hadronic shower prong multiplicity dispersion as a function of (rich) for u-A and u-A 
interactions. 



For charged particle production in hadronic interactions, however, an empirical relation ll33l was found 

D,^ = A + BM. (5) 

The linear dependence of the D^h on (rich) would imply that the shape of the distribution of the 
normalized multiplicity is independent of the hadronic effective mass W. Koba, Nielsen and Olesen have 
shown that at asymptotically high energies, a probability P(nch) to produce an interaction of multiplicity 
rich scaled by the average (rich) is independent of the energy. 

where is the probability distribution independent of energy. This is called Koba-Nielsen-Olesen (KNO) 
scaling. Exact KNO-scaling demands that, in the linear dependence of the dispersion on the average 
multiplicity, Eq. [51 the intercept parameter A is exactly equal to zero. For our event samples, D^h is 
given as a function of {rich) in Fig. |4] with a linear fit superimposed. The values of the fit parameters are 

Dchii^-A) = (1.18 ± 0.17) + (0.20 ± 0.05)(nch) 

Dch{T^-A) = (1.16 ± 0.21) + (0.24 ± 0.08)(nch) 

The fitted values of A and B obtained in this and other experiments are shown in TablefTOl One can 
observe that our data, as well as freon data, aie not compatible with A = 0. Indeed, A = Oh expected 



Table 7: The values of the parameters a and b obtained by fitting the relation rich = « + blnW^ to 
the average charged hadron multiplicity rich- The results of other neutrino experiments are shown for 
comparison. 'Em' stands for nuclear emulsion; (n+) and (n_) are the mean multiplicities of positively 
and negative charged particles. 
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Table 8: (QE-i-RES)-like events in u-A and v-A interactions. 



ng 1 Total 

nch = 39 39 78 

'-'ti rich = 1 56 23 79 

Vf, nch = 132 25 157 

Vf, rich = 1 44 44 



only for a proton or neutron target. A non zero value of A is due to the heavy nuclear targets in nuclear 
emulsions. In such a case one introduces a new variable z' defined as [9] 



rich - a 



(7) 



Here a 



{rich - a) 

(no) and (no) is the extrapolated point where the fitted dispersion line crosses the 



• ]<W<]8 




0.5 1 1.5 2 2.5 

z'=(n . -a)/(<ii,. >-a) 




z=(n^^-a)/(<n^.^>-a) 



Figure 5: The KNO scaling distribution for shower prongs. The superimposed curve represents a fit to pp 
data 1 34]. The data approximately agree with KNO scaling, i.e., the data points at different W'^ intervals 
lie approximately on a single curve. 



Table 9: Kinematic variables for (QE+RES)-like events for and i/^. 
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Table 10: Parameters of the linear fit to the dispersion distributions. The results of other neutrino ex- 
periments are also shown. Note that L>+ and Z)_ are the dispersion for positive and negative charged 
particles of reaction z^-A and I7-A, respectively. 
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average multiplicity axis and found to be -5.82 and -4.83 for u-A and u-A interactions, respectively. 
The KNO-scaling law now reads 

K' = {M -a)PM = ^{z'). (8) 

Clearly, in the asymptotic limit it makes no difference whether one uses z or z'; in fact this expression 
describes an approach to restore KNO scaling. Figure [5] shows the distributions obtained for ^(2;') as a 
function of z' for two intervals of W"^, for v-A and u-A samples; they are compatible with KNO scaling. 
Similar results have been obtained in other experiments |[35l[36l[37ll . however, with different values of a. 

8 Conclusion 

The multiplicity features of (anti-)neutrino-nucleus interactions in emulsion have been investi- 
gated. The results presented in this paper have been obtained with the main objective to aid in tuning 
(anti-)neutrino-nucleus interaction models of the Monte Carlo event generators. To ease the use of these 
results, the numbers are presented in detail in the form of tables. The results can be summarized as 
follows. 

For the first time with nuclear emulsion, the charged particle (shower, grey, black tracks) produc- 
tion at the neutrino interaction vertex of (QE-i-RES)-like events has been studied. It has been found that 
in ( 1.2-1-0. 4-1-0. 2)only the muon track is seen in the final state. These interactions are mainly from the 
reaction u^n n^p and the proton is absorbed by the nucleus without any visible recoil. We also report 
the first study of anti-neutrino induced events in nuclear emulsion for the production of charged parti- 
cles. The average number of shower and heavy prongs in anti-neutrino induced events are measured to 
be (rech(i^-A)) = 2.8 ± 0.1 and (nh(F-A)) = 3.5 ± 0.2, respectively. The fraction of (QE-i-RES)-like 
events in anti-neutrino interactions is found to be (26.3±1.4±3.9)%. 



The dependence of the average multipUcity (rich) on In W"^ for u-A and V-A interactions is com- 
patible with being linear with similar slopes. The dispersion Dch of the multiplicity distribution shows 
a linear dependence on mean multiplicity (rich)- The emulsion data are consistent with the KNO scahng 
as a function of an appropriate multiplicity variable z'. 
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